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Designing appropriate techniques to effectively control the trapping process in complex systems towards 
desirable efficiency is of paramount importance in the study of trapping problem. In this paper, we present 
three different methods guiding trapping process in a scale-free small-world network with a deep trap 
positioned at an initial node. All the proposed approaches dominate the trapping process by varying the 
transition probability of random walks. In the first two techniques, the transition probability is modified by 
an introduced weight parameter and a stochastic parameter, respectively. And the third scheme is a 
combination of the first two approaches, controlled by both parameters synchronously. For aU the three 
control strategies, we derive both analytically and numerically the average trapping time (ATT) as the 
measure of the trapping efficiency, with the obtained explicit expressions being in good agreement with their 
corresponding exact numerical solutions. Our results indicate that the weight parameter changes 
simultaneously the dominating scaling of ATT and its prefactor. Different from the weight parameter, the 
stochastic parameter only modifies the prefactor, keeping the leading scaling unchanged. Finally, compared 
with the first two manners, the third strategy is a fine control, possessing the advantages of the first two ones. 
This work deepens the understanding of controlling trapping process in complex systems. 

Trapping is a paradigmatic dynamical process, which can describe a variety of other dynamical processes 
occurring in diverse complex systems, e.g., lighting harvesting'"^, energy transport' ' and target searching"''. 
Trapping problem constitutes an integral primary problem of random walks, defined as a kind of random 
walks taking place in networks in the presence of a perfect trap, absorbing all particles that visit it'°. A basic 
quantity relevant to the trapping problem is trapping time, also called mean first-passage time (MFPT)""". The 
MFPT from a node i to the trap is the expected time taken by a particle leaving from / to reach the trap for the first 
time. The average trapping time (ATT) is defined as the average of MFPTs over all starting nodes other than the 
trap. The ATT is a primary indicator of trapping, since it measures the efficiency of trapping process. Hitherto, 
evaluating ATT has received considerable attention for diverse networks'' 

Most previous works focused on unraveling the effects of geometrical structure on trapping efficiency. In the 
research area of trapping, another central problem is to control the dynamical process. Recently, controllability of 
complex networks has become a topic of active pursuit and has attracted extensive attention'" Within the 
framework of controllability, several basic issues related to network control have been addressed, including 
optimization of controllability', energy required in control^'', centrality controP\ controlling edge dynamics^'', 
bimodality in control process"", and capacity control"*. In the context of trapping in complex systems, it is 
desirable to design control techniques steering the trapping process towards wanted trapping efficiency. 
Significant efforts have been devoted to this problem in several networks, such as one-dimensional systems^', 
dendrimers', and treelike fractals'"'. However, related control approaches for trapping process in scale-free small- 
world networks remains less understood, in spite of the fact that this network family displays some prominent 
features of real systems'" 

Although real-life networks are often random, the study of artificial networks with well- controlled properties is 
also of high interest*''. In this paper, we focus on the trapping process in a hierarchical small- world scale-free 
network''^, which has received much recent interesf"""^'. We propose three control strategies to steer the trapping 
process with a perfect trap fixed at an initial node. All the three methods control the trapping efficiency by 
modifying the local transition probability. In the first technique, a positive tunable weight parameter is intro- 
duced, which acts as a similar role of energetic funnel in dendrimers''''. The second approach is based on delay- 
diffusion^^'^', where a stochastic parameter is introduced to change the transition probability. The third manner is 
a combination of the first two methods, which is controlled simultaneously by the weight and stochastic 
parameters. 
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For all the three control strategies, we derive analytically the ATT 
and their leading scalings, which are consistent with the numerical 
results. We show that the weight parameter has a strong effect on 
both the leading scaling and the prefactor of the trapping efficiency, 
by changing which the ATT can scale as a superlinear, linear, or 
sublinear function of the system size. The stochastic parameter assoc- 
iated with the delayed diffusion can only modify the coefficient of the 
dominating scaling for ATT, without changing the leading behavior 
of trapping efficiency. Finally, since the third approach combines the 
advantages of the first two methods, it can finely control the trapping 
process. This work opens new avenues to effectively control trapping 
process towards a desirable case with ideal efficiency. 

Results 

Construction and properties of the network. We first introduce the 
construction and some relevant features of the small-world scale-free 
network, which is buUt in an iterative manner. Let F„ (« s 0) denote 
the network after n iterations (generations). Then, the construction 
of F„ can be described as follows. Initially {n = 0), Fq contains two 
nodes connected by an edge. For any « > 0, F„ is obtained from F„ _ i 
by performing the following operation. For each edge in F„ _ i , a new 
node is introduced and linked to both end nodes of the edge. Figure 1 
shows the construction algorithm schematically. 

The particular construction allows determining relevant prop- 
erties of F„. The number of newly introduced nodes at generation / 
(/ s 1) is N, = 3'^^ Thus, the number of nodes N„ in F„ is 

N J 

N„=2+^N,= -(3" + 3). (1) 
i = i ^ 

Let Q.„ represent the set of nodes in F„, and let Q„ be the set of nodes 
created at generation n. It is evident that n„ = Q„_i Ui^n> I^^kI ~ 
= (3" + 3)/2, and |n„| = JV„ =3"^\ where |Q| denotes the number 
of elements in set Q. 

In network F„, the degree of all nodes created at the same iteration 
is identical. Let fc,(n) stand for the degree of node in F„, which fulfils 
the relation 

k,{n+l)=2k,{n). (2) 

Let K„ denote the sum of degrees for all nodes. Then, 

K„=2K„_i+2N„, (3) 




n = 3 

Figure 1 | Illustration for the evolution process of the small-world scale- 
free network. 



which, together with the initial condition Kg = 2, leads to 

K„=2x3". (4) 

Thus, the average degree of F„ is {k)„ = K„/N„ = (4 X 3""')/(3""' + 
1), which is approximatively equal to 4 when the network is very 
large, indicating that the network is sparse. 

The network being studied presents some remarkable properties 
as observed in most real systems"'^'*'''"'. It is scale-free with the degree 
of nodes following a power-law degree distribution P{k) ~ k~", where 
the exponent y = 1 + In 3/ln 2. Its average path length, defined as the 
mean of the shortest distance over aU pairs of nodes, increases loga- 
rithmically with the network size'*'" '". Moreover, in the large network 
size limit, its average clustering coefficient tends to 4/5. Thus, the 
network exhibits the small-world effect'"''. 

Formulation of trapping in the network with a perfect trap at an 
initial hub node. As shown above, the goal of the present paper is to 
introduce three techniques to control the trapping process in the 
scale-free small-world network F„. Since any trapping process is 
dominated by the transition probability from one node to another, 
below we will address three control approaches by modifying the 
transition probability in different ways. In the first method, we 
change the transition probability by introducing a weight for every 
edge, which leads to qualitatively different behavior in the scaling of 
ATT; in the second approach, we modify the transition probability 
according to the birth times of different nodes, which cannot change 
the leading scaling of ATT but can vary its prefactor. The third 
manner is a combination of the former two approaches, which can 
alter the ATT both qualitatively and quantitatively. 

Before introducing the control methods, we first formulate the 
problem of trapping in F„ in the presence of a deep trap located at 
an initial node. We focus on the discrete-time trapping problem 
described by discrete-time random walks. For the convenience of 
description, we label all the nodes in F„ in the following way. The 
initial two nodes in Fq are labeled as 1 and 2, respectively. In each new 
generation, only the new nodes created at this generation are labeled, 
while the labels of all old nodes remain unchanged, i.e., we label new 
nodes in F„ as N„_i + 1, ]V„_i -|-2, ■ ■ ■ , N„. Without loss of gen- 
erality, we assume that the trap is placed at the initial node labeled by 
1. Note that throughout the whole paper, we only consider the near- 
est-neighbor random walks, where the walker is not allowed to stay 
put. At each time step, the probability of the walker moving from its 
current position, say node / (/ s 2), to one of its neighboring nodes j (j 

7^ i) is pij, satisfying 2__^ ^ i Py ~ 1- Moreover, for two different adja- 
cent nodes j and k of node i, the transition probability might be not 
equal to pik- 

(n) 

Let T- stand for the trapping time, that is, the expected time for a 
walker starting from node i to arrive at the trap in F„ for the first time. 

By definition, t'"' obeys the following relation 

Equation (5) characterizes the Markovian property of the random 
walks, which can be recast in matrix notation as 

T = PT-he, (6) 
where P = [py] is a matrix of order N„ — 1, 
T= (^T2"\tI"\ ■ ■ ■ Tn"'^ is an (N„ — 1) -dimensional vector, and 

eisthe(N„ — l)-dimensionalvectorofones,i.e.,e= (1, 1, • • • , 1)^. 
From equation (6), we obtain 
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(I-P)"'e, 



(7) 



where I is the {N„ — 1) X {N„ — 1) identity matrix. 

Using equation (7), the ATT, denoted by (T)„, which is the mean 
of over all non-trap starting nodes distributed uniformly over 
the whole network F„, can be obtained by 



1 



N„ 



(8) 



2 ; = 2 



where Ty is the ijth element of matrix (I— P)"'. The quantity {T)„ is 
significantly important since it is a quantitative indicator measuring 
the trapping efficiency: the less the ATT, the more efficient the trap- 
ping process, and vice versa. 

Equation (8) shows that the problem of determining ATT (T)„ can 
be reduced to finding the sum of all elements of a matrix. Although 
the computation of inverting a matrix puts heavy demands on time 
and memory for large networks, equation (8) can be used to check the 
results for ATT derived by other methods, at least for networks with a 
small number of nodes. In what follows, we will apply equation (8) 
to verify our analytical solutions for the three proposed control 
techniques. 

Qualitatively controlling the trapping efficiency by introducing 
edge weight. The first control method changes the local transition 
probability through introducing a weight for each edge. As will be 
shown, using this method, the ATT can behave superlinealy, linearly, 
or sublinearly, with the system size. 

Numerical solution to ATT in a weighted scale-free small-world net- 
work. The aforementioned binary network F„ can be extended to a 
weighted network, by introducing a weight parameter w (w > 0) in 
the construction algorithm. We use W„ to denote the weighted net- 
work, which is built as follows. For « = 0, Wo is composed of two 
nodes linked by an edge with unit weight. For n s 1, W„ is derived 
from W„ _ 1 in the following way. For each existing edge in W„ _ i , a 
new node is created and attached to both end nodes of the edge by 
two new links. The weight of both new edges is w times that of the 
original edge giving birth to the two new edges. Repeating the iter- 
ative steps, we obtain the weighted version of F„. Note that W„ 
reduces to the binary case F„ when w = I. 

By construction, the total strength S„ of network W„, namely the 
sum of weights of all edges, follows the recursive relation: 

S„=S„_i-|-2wS„_i. (9) 



Considering Sq = 1, we have 



S„ = {lw+\) 



(10) 



Mathematically, the properties of network W„ are reflected in the 
generalized adjacency matrix (weight matrix) W with the element tVy 
specifying the weight of the edge linking nodes / and;'. Let s, denote 
the strength^^^ of node i in W„, which is the sum of weights over 

._^Wij. Notice that 

s, is the fth nonzero entry of the diagonal strength matrix 
S = diag(si, S2, ■■■ , sjv J of network W„. 

After defining the weighted network W„, we now study trapping 
in W„ with a deep trap positioned at an initial node 1. For discrete- 
time random walks in W„, at each time step, the walker moves from 
its current location, say node i, to one of its neighboring nodes j by 
probability = w^/s,. In fact, all the entries py (2 S j s JV„) 
constitute matrix P that is an (]V„ — 1) X (N„ — 1) sub-matrix of 
S^'W with the row and column corresponding to trap being 
removed. Thus, the ATT (T)„ can be computed by using equation (8). 



In addition to equation (8), the exact value of ATT {T)„ for trap- 
ping in W„ can be also determined in terms of the eigenvalues and 
eigenvectors of the Laplacian matrix of network W„, by applying the 
universal framework for random walks in weighted networks^^''. We 
use L to represent the Laplacian matrix of W„, defined as L = S — W. 
Let /li, /,2, ■ ■ ■ , ^N„ be the N„ eigenvalues of L, rearranged as 
0 = /ti </l2 < ■ ■ ■ <In„, and let Hi, ji2, ■ ■ ■ , IJ-m,, be their corres- 
ponding mutually orthogonal eigenvectors of unit length, where 
lij= iJ.j2, ■ ■ ■ , Then, the quantity (T)„ can be repre- 

sented as 



N„ 



N„ 



- k = l 
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We have performed extensive numerical computation for (T)„ using 
both equations (8) and (11), the results of which are consistent with 
each other. In the following text, we will utilize these numerical 
results to check our analytical expression for (T)„. 

Although one can employ equations (8) or (11) to obtain numer- 
ically but exactly the ATT (T)„, either inverting a matrix or deter- 
mining the full spectrum of a matrix is prohibitively time and 
memory consuming, making it intractable to employ equations (8) 
or (11) to compute (T)„ for a large network. In particular, by making 
use of equations (8) or (11), it is hard and even impossible to obtain 
asymptotic scaling for the ATT (T)„. Hence, it is of utmost import- 
ance to develop a computationally cheaper method for computing 
(T)„. Fortunately, the peculiar construction and structure of the 
studied network allow to analytically treat (T)„ as well as the ATT 
for the other two trapping problems, obtaining closed-form solutions. 

Analytical solution to ATT. For the purpose of studying ATT ana- 
lytically, let's begin with examining random walks in W„ + 1 . By con- 
struction, upon growth of the network from generation n to next 
generation n + I, the degree of node / doubles, that is, it grows from 
ki{n) to ki{n + 1) = 2A:,(n). These 2fc,(«) neighbors of ; can be 
classified into two categories: one half are the original neighbors 
belonging to Q.„, and the other half are newly introduced neighbors 
in n„ + 1 . Let Z represent the MFPT for going from node i to any of its 
kj{n) old neighbors; and let Xbe the MFPT for going from any of the 
ki(n) new neighbors of / to one of its /c,(n) old neighbors. Then we can 
establish the following underlying backward equations: 

1 w 
— + y^(l+X), 
l + w l + w (^2) 

X- 2 + 2^^+^)' 

which lead to 2 = (4w + 2)l{w + 2). Thus, upon the evolution of the 
weighted network from generation n to generation « -I- 1, the trap- 
ping time for an arbitrary node ; increases by a factor of (4w -I- 2)1 {w 
+ 2), that is, 

r(" + l) 



4W-\-2 ^(„) 



(13) 



which will be helpful for the following evaluation of the exact formula 
for<I)„. 

Having obtained the evolution rule of trapping time for any node 
in the weighted network, we continue to determine the ATT in W„. 
For this purpose, we define the following intermediary quantities for 
1 <g<n: 

ien. 



and 



(„) 



(15) 
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Then, (T)„ is given by 



iT)„- 



' ^ «,tot' 



(16) 



N„-l 

which reduces the problem of evaluating (T)„ to determining r^"o,. 
By definition, 

„(«) _^(«) ^(n) _4w + 2 („_i) -(„) 
n.tot ~ n- l,tot "T «,tot ^ I T n — l.tot "T K.tot' 

w-\- Z 

where equation (13) is used. Equation (17) indicates that, in order to 
obtain T^"o(, we should first deduce f^"!)!- 

By construction, for each new node x in f2„, it has only two neigh- 
bors u and V, both belonging to n„_i. The ATTs of these three nodes 
obey 

T(«) = i + ir(«) + irw. (18) 

Summing equation (18) over all nodes in Q„, we obtain 



f(") =|o 1+ V 

^ n,tot 1 I ^ / J 

ien„_i 



fc,(«-i)x^r^ 



(«) 



based on which we can further derive 

r("+i) - In I ^ 



h{n)y.-^T\ 



(n+l) 



(19) 



(20) 



Equation (20) minus equation (19) times 2(4vi/ + T)l(w + 2) and 
recalling equations (2) and (13), we obtain the following recursive 
relation: 



^(« + i) _ 

«+ l.tot 



8w + 4 - 

T 

w + 2 

8^ + 4 



.(„) 



w + 2 



(n+l) 



(21) 



-lo I 8w+4 4w+2 („) 

Using |fi„| =3""' and f\]l^ = 3{w+ 1) /{w + 2), equation (21) can 
be resolved inductively to yield 

3" " / 6 

(llw-l-4) 4 I +10^-4 



r(") — 

^ «,tot ~ 



2w 



w + 2 



22 



Substituting equation (22) into equation (17) leads to 

(„) _ 4w+2 („_i) 



w + 2 

^n — 3 
21^ 



;ilw + 4) 4 



w + 2 



-l-lOw-4 



(23) 



Considering ,2,, = 1 -|-3(2m'-|- 1)/ (1^ + 2), equation (23) is solved to 
yield 

3«-i 



T<"'o.= ^(2x3" + « + 4) 



24 



and 



r(«) 



36^(^-4) 



[(llM' + 4)(tv-4) 



6 + l2w 
w + 2 



(25 



+ (451^^-72^)1^^^) +3"(16-32w-20m^)] 



for w = 4 and w 4, respectively. 

Plugging the results in equations (24) and (25) into equation (16), 
we arrive at the closed-from expressions for (T)„ given by 



and 



_ 2x3" + n + 4 
(J)„- 3„^3 X3 



1 /4w + 2 

, r[(llw + 4)(tv-4)3" — 

18w(iv-4)(3" + 3) ' \w + 2 



+ (45^^-721^) I I +3"(l6-32w-20w^^ 

w + 2 



(26) 



(27) 



for w = A and w 4, respectively. When w = 1, W„ reduces to the 
unweighted network F„, and equation (27) is consistent with the 
result previously obtained in Ref 57. 

We have checked our analytical formulas against numerical results 
obtained from equations (8) and (11). For different values of w and n, 
the analytical results obtained from equations (26) and (27) comple- 
tely agree with those numerical ones; see Fig. 2. This agreement 
serves as an independent test of our theoretical formulas. 

Equations (26) and (27) show that for trapping in W„, the exact 
expression of ATT to an initial hub node is dependent on the weight 
parameter w. We now proceed to show how the leading behavior of 
ATT scales with the number of nodes in network W„. 

Leading scalings of ATT. Recalling N„ = (3" + 3)/2, equation (26) 
shows that (T)„ ~ N„ for w = 4. When w 4, the first term on the 
right-hand side (rhs) of equation (27) dominates, implying that 
(r)„~N>(''"'+2)/("'+2). Therefore, in the whole range of w > 0, 
the leading scaling of (T)„ can be unified as 



(28) 



4W-I-2 

where the exponent ii{w) ^log, — is an increasing function of 

parameter w and can be less than, equal to, or greater than 1. 
Concretely, if 0 < w < 4, rjiw) < 1, implying that (T)„ scales sub- 
linearly with N„; if w = 4, r]{w) = 1, implying that (T)„ behaves 
linearly with N„; and if w > 4, ii{w) > 1, implying that (T)„ varies 
superlinearly with N„. 

Equation (28) indicates that the parameter w has a significant 
impact on the trapping efficiency, by adjusting which the ATT can 
display rich behavior revealed by exponent rjiw). When w grows 
from zero to infinite, iliw) increases monotonously from 0 to 
log3 4. Thus, one can adjust the parameter w to modify the transition 
probability and qualitatively control the trapping process. 



10000 -. 




Figure 2 | ( 7)„ as a function of n for three different values of w. The filled 
symbols are the data coming from numerical results obtained by direct 
calculation from equations (8) and (11); while the empty symbols 
correspond to the exact analytical values given by equations (26) or (27). 
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Quantitatively modifying the trapping efficiency by delayed ran- 
dom walks. We now present a control technique to quantitatively 
alter the trapping efficiency. It is often the case, when a network 
assembly is finished, its new topology does not come into play in 
dynamical processes right now. That is, there is delay for the new 
structure to influence dynamical processes taking place on the 
network. This phenomenon inspires us to propose the second 
control strategy that includes the delay of operation for new 
structure. 

Definition of delayed random walks. By construction, the structure of 
F„_i is embedded in F„. Then, we can define the following random 
walks with delay, called delayed random walks, performed on the 
binary network F„. During the process of delayed random walks, a 
walker starting from an old node in F„, which is already existent in 
F„_i, is allowed to perform isotropic nearest-neighbor random 
walks in either F„_i or F„, with respective probabilities p and 1 — p 
(0 S p < 1). That is to say, for delayed random walks in F„, the 
transition probability from an old node to one of its adjacency nodes 
j is given by 



(29) 



If the current state of the walker is at a new node, then it performs 
uniform nearest- neighbor random walks in F„. 

There are two limiting cases for the above-defined delayed random 
walks in F„. Whenp = 0, it reduces to unbiased random walks in F„; 
when p = 1, it is exactly unbiased random walks in F„_i. As 
expected, the probability parameter p governs the process of delayed 
random walks in F„ . Below we will show that by varyingp, significant 
modification occurs in the prefactor of ATT, with its dominating 
behavior unchanged. 

Closed-form expression for ATT. We here concentrate on delayed 
random walks in F„ with a perfect trap at an initial node. In this case, 
we use F,- to represent the trapping time for a walker starting from 
node i to first arrive at the trap in F„, and use to denote the ATT. 
In order to determine we define two quantities for ^ < n: 

^It = i:,en, ^""^ ^i"o.= Then, the ATT {F)„ is 

given by 



N„-l 



Next, we determine P, 



(") 

An) 



(30) 



(«) 



Before evaluating F„ tm, we first derive the evolution rule of 
Analogous to random walks in W„ + i, the quantities Z and X for 
delayed random walks follow relations: 



Z=p+{l-p) 
1 1 



x= 



2 2 



11/ 

(1+2), 



(31) 



which lead to Z = 2(3 — p)/(3 -I- p). Thus, we have 

2(3-p) 



3+p 



(32) 



where T^"^ is trapping time of node ; for random walks in 
corresponding to the particular case of w = 1. 



Using equation (32), F^"ot can be computed as 



p(") _ pi") 



(«) _ 2(3-p) ^(„_i) _^-(„) 



pi"; _ 

-■^^n.tot ~ 



3+p 



(33) 



Substituting w = 1 into equation (25) gives r^"_i tit equation (33) 



t: 



{„-!) 



= — (5x6"-'-h3x2"-'+4x3"-') 
12 ^ ^ 



(34) 



Thus, to obtain P^^iov we only need to determine P^^iot- 

By using a similar process as that for random walks in ^ 
obtain 



p(") - lO I 4_ 
-^n.tot ~ ]"nl "T / J 



in) 



fcj(«-l) X -P\ 



/c,(„-i)x^r("-'' 



and 



hin + l) 



3-hp ' 



(35) 



3+p ' 



(36) 



Equation (36) minus equation (35) times 4 and considering the 



relations T^"' —2T-" andA:,(«) = 2kj{n — 1), we obtain the recurs- 
ive relation: 



r(«-i 



^n+l,tot ^^n.tot" 



— i^J7+l 



-4 Q„ 



2(3 -p) 
3 -hp 



.(«) 



2(3 -p) ^{„) 
3+p 



(37) 



Combining equation (23) for w = 1, |n„|=3" ^ and the initial 
condition p[^^^^ = 6/{p + 3),v/e can resolve equation (37) to obtain 

on — 2 



pi") 



2(3 +p) 



[5(3-p) x2" + 10p-h6]. 



(38) 



Plugging equations (34) and (38) into equation (33) leads to the 



explicit expression for pl^tg^: 



tU'^x2"-^+ y~^i 6"-'+3"-', 
3 -hp 2(3 -hp) 



(39) 



which, together with equation (30), produces the exact formula for 
the ATT {P)„ given by 



6(3" -h 3) 



'"^x2"-^+'i^~^j6"-'+3"-' 



3+P 



2(3 +p) 



(40) 



We have confirmed the analytical solution in equation (40) by com- 
paring it with equation (8). The analytical and numerical results are 
completely consistent with each other, indicating that the predicted 
formula in equation (40) is valid. In addition, for the two special cases 
of p = 0 andp = 1, equation (40) agrees with the values of (T)„ and 
for random walks in W„ corresponding to tv = 1, also indi- 
cating the validity of equation (40). 

Leading behavior ofATT and its prefactor. Considering N„ = (3" -h 3)/2, 
the expression for {P)„ in equation (40) can be represented in terms of 
the system size N„ in the following form: 



1 



-h 



12N, 
3-p 



5(3 -p) 
4(3 +p) 



(2N„-3)'+'°8^' 



,(2N„-3)'°832+^N„-l 
4(3 +p)^ ^ 3 



(41) 
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which provides the exact dependence of (P)„ on the network size JV„ and 
parameter p. For a large system, i.e., N„ we have the following 

expression for the dominating term of {P)„: 

5(3 -p) 



24(3 +p) 



(2N„ 



(42) 



Equation (42) indicates that in the whole range of 0 £ p £ 1, the ATT 
increases subUnearly with the system size N„, with the exponent log3 2 
independent of parameter p. Although parameter p has little influence 
on the dominating behavior of ATT, it modifies significantly the pre- 
factor [5(3 — p)]/[24(3 + p)] of the dominating term of ATT, which is a 
descending function of p. When p grows from 0 to 1, the prefactor drops 
from 5/24 to 5/48, suggesting that the inclusion of delayed random walks 
can improve the trapping efficiency in a significant way. For the two 
extreme cases of p = 1 and p = 0, the ATT for the former with maximal 
delay is only one half of that for the latter without any delay. 

Unifying control technique. Thus far, we have described two 
different approaches to control the efficiency for trapping in the 
network concerned. The first method can change the leading term 
of ATT and its coefficient simultaneously, while the second 
technique only modifies the prefactor but keeps the leading term 
unchanged. In the sequel, we will present a unifying control 
scheme incorporating the two previous methods. 

Explicit formula for ATT. The new technique is in fact a delayed 
random walk in the weighted network W„. For this case of random 
walks, if the walker is currently at an old node belonging to n„_i, it 
can perform random walks either in W„_i with probability p or in 
W„ with complementary probability 1 — p; if the walker is currently 
at a new node belonging to fl„, it performs random walks in W„. 
Then the quantities Z and X obey the relations 



Z = p + (l-p)[^ + J^(l+X) 



(1+2), 



(43) 



which give Z = 2(1 + 2w — wp)/{2 + w + wp). 

For delayed random walks in W„ with a trap fixed at an initial 

node, we use P,-"' to denote the trapping time for node /, and use (P)„ 
to represent the ATT to the trap. Then, from equation (43), we have 

(44) 



p(« + i)^ 2(l + 2w-wp) ^(„) 



2 + w + wp 



In order to determine {P)„, ioT g £ n, we define P^,tot = ^ Pj 
efore. Then, 



N„-l 



45 



where P^"ot can be rewritten as 



p(«) _ p(") I p(") 

-'^f7,tot -'^/j — l,tot n.tot 



2(l + 2w-wp) („_i) („) 
2 + w + wp '"-1.""+^"."" 



(46) 



The quantity P^"ot can be evaluated as follows. First, it is easy to 
derive 



p<"' =ln 1+ V 



' 2 + w + wp 



(47) 



and 



p("+i) _ lo I ^ 



fc,(»)x ^+'"-^^ r,'"' 

2 + w + wp 



, (48) 



both of which yield the relation 



(„ + i) _ 2(4w + 2) („) 2(4w + 2) 

W-r z w -\- 

2(l+2w-wp) -(„) 



(49) 



2 + w + vvp 



Using pj'f'g, = (3 + 3w)/(2 + M'+pvi'), equation (49) is solved to 
obtain 



2w{2w+l)[{p+l)w + 2] 



(w + 2)(llw + 4)[(2-p)w+l] 



4^ + 2^ (50) 



w + 2 



+ 2(2^+1) [w(llM' + 4)p + 5w^ + 8w-4] }. 
Substituting equations (24), (25) and (50) into equation (46) yields 



p(") _ 
-^n.tot ~ 



and 



4p + 6 



1 



[(9-4p)(2x3" + n)+4(9+p)] 



(51) 



d(") ^ 

54h'[(p+1)w + 2; 

f3[(2-p)w+l](llM/ + 4)(w + 2) (\2w + 6 



+ 



2(2w+l) V w + 2 

27vf[(2-p)w+l](5vf-8)(w + 2) /4w + 2 



(2w-8)(2w+l) 



w + 2 



(52) 



3" 



H 6(7p-5)w^-12(4p + 9)M'^ 

-24(2p + 3)w + 48]}, 

for w = 4 and w 4, respectively. 

Then, according to equation (45), the explicit expression for ATT 

is 



2x3" 



(4p + 6)(3« + l) 



[(9-4p)(2x3" + n)+4(9+p)] (53) 



and 



54w[(p+l)w + 2](3« + l) 
3[(2-p)w+l](llM' + 4)(w + 2) /I2w + 6y 

2(2w+l) V 1^+2 / 

27w[(2-p)w+l](5w-8)(w + 2) /4w + 2' 



+ 



(2h'-8)(2w+1) 



w + 2 



(54) 



+ 6(7p-5)w -12(4p + 9)tv^ 

-24(2p + 3)w + 48]}, 



for w = 4 and w 7^ 4, respectively. 

For the particular case w = 1, equation (54) reduces to equation (40); 
and for the special case p = 0, equations (53) and (54) reduce to 
equations (26) and (27), respectively. Thus, for these two limiting 
cases, equations (53) and (54) are correct. Furthermore, we have also 
verified the exact expressions given by equations (53) and (54) 
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against direct computation from equation (8). For various w andp, 
the analytical and numerical results agree with each other. 

Dominating scaling for ATT and its prefactor. For a very large net- 
work, according to equations (53) and (54), the leading term of (P)„ 
can be expressed in terms of network size N„ in a single formula as 

_(2-p)m/+1(11w + 4)(w + 2) 



(l+p)w + 2 



18w 



2N„ 



55 



From equation (55), we can see that the weight parameter w deter- 
mines the dominating scaling of {P)„, and that parameters w and p 
together determine the prefactor. Therefore, in order to obtain 
needed trapping efficiency, one can first adjust parameter w to attain 
the desired dominating scaling, then continue to change parameter p 
until the prefactor is desired. 

Discussion. We have proposed three techniques controlling the 
efficiency for trapping in a small-world scale-free network with an 
immobile trap located at an initial node. The first method can 
qualitatively change the leading scaling of trapping efficiency, 
which is realized by introducing a parameter governing the weight 
of each edge. By tuning the weight parameter, the ATT can scale 
superlinearly, linearly, or sublinearly with the system size. The 
second approach is based on delayed random walks, controlled by 
a stochastic parameter that has little influence on the leading scaling 
of ATT but can significantly affect its prefactor. The third method is a 
combination of the firth two approaches, it can thus vary both the 
leading scaling of ATT and its prefactor. For the three trapping 
problems, we have studied both numerically and analytically the 
ATT, the results of which are in good agreement with each other. 
We expect that similar techniques can also be applied to control other 
dynamics performed in complex networked systems, such as 
synchronization^**, navigation""', and search'". 



1. Bar-Haim, A. & Klafter, J. Geometric versus energetic competition in light 
harvesting by dendrimers. /. Phys. Chem. B 102, 1662-1664 (1998). 

2. Bar-Haim, A. & Klafter, J. Dendrimers as light harvesting antennae. /. Lumin. 
76-77, 197-200 (1998). 

3. Bar-Haim, A., PQafter, J. & Kopelman, R. Dendrimers as controlled artificial 
energy antennae. /. Am. Chem. Soc. 119, 6197-6198 (1997). 

4. Bentz, J. L., Hosseini, F. N. & Kozak, J. J. Influence of geometry on light harvesting 
in dendrimeric systems. Chem. Phys. Lett. 370, 319-326 (2003). 

5. Bentz, J. L. & Kozak, J. J. Influence of geometry on light harvesting in dendrimeric 
systems. 11. nfh-nearest neighbor effects and the onset of percolation. /. Lumin. 
121, 62-74 (2006). 

6. Sokolov, I. M., Mai, J. & Blumen, A. Paradoxal diffusion in chemical space for 
nearest-neighbor walks over polymer chains. Phys. Rev. Lett. 79, 857-860 (1997). 

7. Blumen, A. & Zumofen, G. Energy transfer as a random walk on regular lattices. 
;. Chem. Phys. 75, 892-907 (1981). 

8. Jasch, F. & Blumen, A. Target problem on small-world networks. Phys. Rev. E 63, 
041108 (2001). 

9. Shlesinger, M. F. Mathematical physics: Search research. Nature 443, 281-282 
(2006). 

10. MontroU, E. W. Random walks on lattices. III. Calculation of first passage times 
with application to exciton trapping on photosynthetic units. /. Math. Phys. 10, 
753-765 (1969). 

1 1 . Redner, S. A guide to first-passage processes (Cambridge University Press, Oxford, 
2001). 

12. Noh, J. D. & Rieger, H. Random walks on complex networks. Phys. Rev. Lett. 92, 
118701 (2004). 

13. Condamin, S., Benichou, O. & Moreau, M. First-passage times for random walks 
in bounded domains. Phys. Rev. Lett. 95, 260601 (2005). 

14. Condamin, S., Benichou, O. & Klafter, J. First-passage time distributions for 
subdiffusion in confined geometry. Phys. Rev. Lett. 98, 250602 (2007). 

15. Condamin, S., Benichou, O. & Moreau, M. Random walks and brownian motion: 
A method of computation for first-passage times and related quantities in 
confined geometries. Phys. Rev. E 75, 021111 (2007). 

16. Condamin, S., Benichou, O., Tejedor, V., Voituriez, R. & Klafter, J. First-passage 
times in complex scale-invariant media. Nature 450, 77-80 (2007). 

1 7. Garza-Lopez, R. A. & Kozak, J. J. Invariance relations for random walks on square- 
planar lattices. Chem. Phys. Lett. 406, 38-43 (2005). 



18. Garza-Lopez, R. A., Linares, A., Yoo, A., Evans, G. & Kozak, J. J. Invariance 
relations for random walks on simple cubic lattices. Chem. Phys. Lett. 421, 
287-294 (2006). 

19. Kozak, J. J. & Balakrishnan, V. Analytic expression for the mean time to 
absorption for a random walker on the Sierpinski gasket. Phys. Rev. E 65, 021 105 
(2002). 

20. Bentz, J. L., Turner, J. W. & Kozak, J. J. Analytic expression for the mean time to 
absorption for a random walker on the Sierpinski gasket. II. The eigenvalue 
spectrum. Phys. Rev. E 82, 011137 (2010). 

21. Wu, B., Lin, Y., Zhang, Z. Z. 8c Chen, G. R. Trapping in dendrimers and regular 
hyperbranched polymers. /. Chem. Phys. 137, 044903 (2012). 

22. Lin, Y. & Zhang, Z. Z. Influence of trap location on the efficiency of trapping in 
dendrimers and regular hyperbranched polymers. /. Chem. Phys. 138, 094905 
(2013). 

23. Lin, Y., Wu, B. & Zhang, Z. Z. Determining mean first-passage time on a class of 
treelike regular fractals. Phys. Rev E 82, 031140 (2010). 

24. Zhang, Z. Z., Guan, J. H., Xie, W. L., Qi, Y. & Zhou, S. G. Random walks on the 
Apollonian network with a single trap. EPL 86, 10006 (2009). 

25. Zhang, Z. Z., Yang, Y. H. & Gao, S. Y. Role of fractal dimension in random walks 
on scale-free networks. Eur. Phys. /. B 84, 331-338 (2011). 

26. Zhang, Z. Z., Xie, W. L., Zhou, S. G., Li, M. & Guan, J. H. Distinct scalings for mean 
first-passage time of random walks on scale-free networks with the same degree 
sequence. Phys. Rev.ESO, 061111 (2009). 

27. Lin, Y., Julaiti, A. & Zhang, Z. Z. Mean first-passage time for random walks in 
general graphs with a deep trap. /. Chem. Phys. 137, 124104 (2012). 

28. Lin, Y. & Zhang, Z. Z. Mean first-passage time for maximal-entropy random walks 
in complex networks. Sci. Rep. 4, 5365 (2014). 

29. Peng, X. Sf Zhang, Z. Z. Maximal entropy random walk improves efficiency of 
trapping in dendrimers. /. Chem. Phys. 140, 234104 (2014). 

30. Liu, Y.-Y., Slotine, J.-J. St Barabasi, A.-L. Controllability of complex networks. 
Nature 473, 167-173 (2011). 

31. Yuan, Z. Z., Zhao, C, Di, Z. R., Wang, W.-X. & Lai, Y.-C. Exact controllability of 
complex networks. Nat. Commun. 4, 2447 (2013). 

32. Liu, Y.-Y., Slotine, J.-J. & Barabasi, A.-L. Observability of complex systems. Proc. 
Natl. Acad. Sci. 110, 2460-2465 (2013). 

33. Wang, W.-X., Ni, X., Lai, Y.-C. & Grebogi, C. Optimizing controllability of 
complex networks by minimum structural perturbations. Phys. Rev. E 85, 0261 15 
(2012). 

34. Yan, G., Ren, J., Lai, Y.-C, Lai, C.-H. & Li, B. W. Controlling complex networks: 
How much energy is needed? Phys. Rev. Lett 108, 218703 (2012). 

35. Liu, Y.-Y., Slotine, J.-J. 8{ Barabasi, A.-L. Control centrality and hierarchical 
structure in complex networks. PLoS ONE 7, e44459 (2012). 

36. Nepusz, T. & Vicsek, T. Controlling edge dynamics in complex networks. Nat. 
Phys. 8, 568-573 (2012). 

37. Jia, T. et al. Emergence of bimodality in controlling complex networks. Nat. 
Commun. 4, 2002 (2013). 

38. Jia, T. & Barabasi, A.-L. Control capacity and a random sampling method in 
exploring controllability of complex networks. 5c/. Rep. 3, 2354 (2013). 

39. Bar-Haim, A. & Klafter, J. On mean residence and first passage times in finite one- 
dimensional systems. /. Chem. Phys. 109, 5187-5193 (1998). 

40. Wu, B. 8f Zhang, Z. Z. Controlling the efficiency of trapping in treelike fractals. 
/. Chem. Phys. 139, 024106 (2013). 

41. Albert, R. & Barabasi, A.-L. Statistical mechanics of complex networks. Rev. Mod. 
Phys. 74, 47-97 (2002). 

42. Newman, M. E. J. The structure and function of complex networks. SIAM Rev. 45, 
167-256 (2003). 

43. Newman, M. E. J. Networks: an introduction (Cambridge University Press, 
Oxford, 2010). 

44. Boettcher, S., Singh, V. & Ziff, R. M. Ordinary percolation with discontinuous 
transitions. Nat. Commun. 3, 787 (2012). 

45. Dorogovtsev, S. N., Goltsev, A. V. 8f Mendes, J. F. F. Pseudofractal scale-free web. 
Phys. Rev. E 65, 066122 (2002). 

46. Bollt, E. M. & ben- Avraham, D. What is special about diffusion on scale-free nets? 
New I Phys. 7, 26 (2005). 

47. Rozenfeld, H. D., Havlin, S. & ben- Avraham, D. Fractal and transfractal recursive 
scale-free nets. New J. Phys. 9, 175 (2007). 

48. Rozenfeld, H. D. & Ben-Avraham, D. Percolation in hierarchical scale-free nets. 
Phys. Rev E 75, 061102 (2007). 

49. Zhang, Z. Z., Zhou, S. G. & Chen, L. C. Evolving pseudofractal networks. Eur. 
Phys. I B 58, 337-344 (2007). 

50. Zhang, Z. Z., Liu, H. X., Wu, B. & Zhou, S. G. Enumeration of spanning trees in a 
pseudofractal scale-free web. EPL 90, 68002 (2010). 

51. Hasegawa, T. & Nemoto, K. Hierarchical scale-free network is fragile against 
random failure. Phys. Rev. E 88, 062807 (2013). 

52. Friesecke, G. Convergence to equilibrium for delay-diffusion equations with small 
delay. /. Dyn. Diff. Equal 5, 89-103 (1993). 

53. Hale, J. 8c Lunel, S. V. Effects of small delays on stability and control (Birkhauser 
Basel, 2001). 

54. Watts, D. J. & Strogatz, S. H. Collective dynamics of small- worldnetworks. Nature 
393, 440-442 (1998). 

55. Barrat, A., Barthelemy, M., Pastor-Satorras, R. & Vespignani, A. The architecture 
of complex weighted networks. Proc. Natl. Acad. Sci. 101, 3747-3752 (2004). 



SCIENTIFIC REPORTS | 4:6274 | DOI: 1 0. 1 038/srep06274 



7 



56. Lin, Y. & Zhang, Z. Z. Random walks in weighted networks with a perfect trap: An 
application of Laplacian spectra. Phys. Rev. E 87, 062140 (2013). 

57. Zhang, Z. Z., Qi, Y., Zhou, S. G., Xie, W. L. & Guan, J. H. Exact solution for mean 
first-passage time on a pseudofractal scale-free web. Phys. Rev. E 79, 021127 
(2009). 

58. Barahona, M. & Pecora, L. M. Synchronization in small-world systems. Phys. Rev. 
Lett. 89, 054101 (2002). 

59. Kleinberg, J. M. Navigation in a small world. Nature 406, 845-845 (2000). 

60. Benichou, O., Loverdo, C, Moreau, M. & Voituriez, R. Intermittent search 
strategies. Rev. Mod. Phys. 83, 81-129 (2011). 



Acknowledgments 



The authors thank Yihang Yang for his assistance in preparing this manuscript. This work 
was supported by the National Natural Science Foundation of China under Grants No. 
11275049 and the National Basic Research Program of China under Grant No. 
2010CB731401. 



Author contributions 

Y.L. and Z.Z.Z. designed the research, performed the research, and wrote the manuscript. 

Additional information 

Competing financial interests: The authors declare no competing financial interests. 

How to cite this article: Lin, Y. Sf Zhang, Z. Controlling the efficiency of trapping in a 
scale-free small-world network. Sci. Rep. 4, 6274; DOI:10.1038/srep06274 (2014). 



This work is licensed under a Creative Commons Attribution-NonCommercial- 
NoDerivs 4.0 Intemational License. The images or other third party material in 
this article are included in the article's Creative Commons license, unless indicated 
otherwise in the credit line; if the material is not included under the Creative 
Commons license, users will need to obtain permission from the license holder 
in order to reproduce the material. To view a copy of this license, visit http:// 
creativecommons.org/licenses/by-nc-nd/4.0/ 



SCIENTIFIC REPORTS | 4:6274 | DOI: 1 0. 1 038/srep06274 



8 



